Optomechanics of liquid crystals for dynamical optical response of photonic structures 
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We show that the mechanical effect of light on the orientational ordering of the crystalline axis of 
a mesophase can be used to control the dynamics of the optical response of liquid crystal infiltrated 
photonic structures. The demonstration is made using a one-dimensional periodic structure whose 
periodicity is broken by the presence of a nematic liquid crystal defect layer. In this study we 
report on output light polarization and/or intensity dynamics that depends on the initial molecular 
ordering and incident light wavelength and intensity. 



PACS numbers: 42.70.Qs,42.70.Df 

I. INTRODUCTION 

Two decades ago, Yablonovitch [1] and John [2] inde- 
pendently introduced the concept of photonic crystals, 
i.e. dielectric systems with spatial periodicity of the re- 
fractive index, which leads to the formation of photonic 
bandgaps. The latter correspond to wavevector ranges 
(hence implying direction and frequency), for which light 
cannot propagate. Indeed, this feature allows to control 
the flow of light [3] provided an appropriate structured 
material properties. To do so, the light itself offers prac- 
tical solutions towards high-resolution three-dimensional 
optical micro- and nano- fabrication techniques [4]. Al- 
though being a necessary first step, structured material 
is however not enough when flexible optical data pro- 
cessing is envisaged, which requires reconfigurable pho- 
tonic crystals. To do so, linear or nonlinear refractive 
index changes have appeared as a straightforward solu- 
tion that can be supplied by external fields of different 
nature (e.g., thermal, electrical, or optical). Obviously, 
the use of light itself is appealing since it naturally bene- 
fits from the strong spatial confinement of the field inside 
the structure, thereby enhancing the optical response and 
paving the way towards all-optical photonic circuitry [5] . 

Among various nonlinear optical materials, the at- 
tractiveness of liquid crystals result from the genuine 
combination of its 'liquid' and 'crystalline' features, 
which allow easy integration into photonic crystal micro- 
architectures and extreme sensitivity to external fields 
at the same time. Such a possibility has been thor- 
oughly studied since the pioneering results of, on the one 
hand, Busch and John, who proposed to use the orien- 
tational ordering of the optical axis of liquid crystals in 
order to tune the properties of photonic bandgaps [6], 
and, on the other hand, Yoshino and coworkers, who 
proposed to benefit from the successive appearance of 
distinct mcsophascs in thcrmotropic liquid crystals when 
temperature is changed [7]. However, neither thermal 
nor electrical tuning of liquid crystals benefits from the 
very nature of a photonic crystal, contrary to the use of 



light itself to trigger refractive index changes. 

The large orientational optical nonlincaritics of liquid 
crystals [8, 9] can be fruitfully exploited to control the op- 
tical response of a photonic structure via the light that 
propagates inside it. Such an issue was addressed in pre- 
vious works in the simplest case of a one-dimensional 
dielectric structure in which is embedded a pure or dye- 
doped nematic liquid crystal layer with initial uniform 
alignment [10-13]. These works, however, have only re- 
ported on static or transient light-induced changes. More 
recently, rotational dynamics of the director (i.e., the 
unit vector n that represents the local average molecular 
orientation of a nematic) has also been experimentally 
demonstrated [14]. 

In this study, we theoretically demonstrate that op- 
tomechanical effects at mesoscale can be used to control 
output light polarization and/or intensity dynamics of 
liquid crystal infiltrated photonic structures. First we 
introduce the system in section II A and the choice of 
two representative light-matter interaction geometries is 
discussed in section II B. Then, the optical threshold 
behavior is addressed in section III. The all-optical dy- 
namical response of the photonic structures is studied in 
section IV and section V concludes the paper. 

II. BACKGROUND 
A. Definitions 

We choose a one-dimensional periodic structure made 
of alternating layers of SiC>2 and TiC>2 with thicknesses 
103 and 64 nm, respectively, which exhibits a bandgap 
in the visible range between 500 and 720 nm. In all our 
simulations, a uniformly aligned nematic liquid crystal 
defect layer with thickness L — 5 /im is located in the 
central part of the structure with five SiC>2 /Ti02 building- 
blocks on each side (see Fig. 1). 

We consider two different kind of initial uniform align- 
ment n for the director. Namely, the initial alignment is 
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either perpendicular (homeotropic anchoring, 'H') or par- 
allel (planar anchoring, 'P') to the nematic slab. Hence, 
no = e 2 or no = e x , where (e x ,e y ,e z ) is the Cartesian 
coordinate system, as illustrated in Fig. 1(a) and 1(b), re- 
spectively. We consider a nematic material with typical 
refractive indices n± = 1.5 and n\\ = 1.7, where symbols 
(_L, ||) refer to directions perpendicular and parallel to n, 
respectively, and Frank elastic constants are chosen such 
that K x /K z = 2/3 and K 2 /K 3 = 1/2. 

The optical properties of these structures are calcu- 
lated in the plane wave approximation, therefore all vari- 
ables depend on coordinate z and time t only. The 
light propagation problem is solved by using the Bcrre- 
man 4x4 matrix approach [15] and taking into ac- 
count that, inside the nematic, the light is coupled to 
the Euler-Lagrange equations that govern the dynam- 
ics of the director [16]. We also introduce the reduced 
spatial coordinate £ = z/L and time r = £/th,p where 
we defined the geometry-dependent characteristic times 
r H = 7l L 2 /(7r 2 A' 3 ) [17, 18] and r P = 7l L 2 /(^ 2 A 2 ) [19], 
71 being the rotational viscosity. 

The director is represented by the usual spherical an- 
gles and $ following n = (sin cos $,sin sin $,cos 0) 
but the different boundary conditions at £ = and £ = 1 
for the H and P geometries impose to distinguish these 
two cases. In the homeotropic case, n(0,r) = n(l,r) = 
e 2 and the director representation is 

nn = (sin ©h cos $h, sin 0h sin <I>h, cos 0h) , (1) 

with 

e H &r) = x;e<rV)sk(™rO> (2) 

m>l 



^,r)^(r) + ^Hr) ^^^ , (3) 

m>l k *> ^' 

m integer, whereas in the planar case the director lies in 
the plane (x, y) with n(0,r) = n(l,r) = e x , hence 

n P = (cos$ P ,sin$p,0) , (4) 

with 

*p (£, r) = *P° ( T ) sin ( m O • (5) 

m>l 

In practice, we retain a large enough number of po- 
lar and azimuthal modes [i.e., sums that appear in 
Eqs. (2,3,5) arc truncated] in order to ensure accurate 
results, which is done following previous works devoted 
to nematic slab alone in the case H [17, 18] and P [19]. 
We introduce the (adiabatic) total phase delay A between 
the extraordinary (e) and ordinary (o) waves due to the 
nematic layer (which turns out to be a measure of the 
reorientation associated with the polar degree od free- 
dom of the director). Again, distinction must be made 
between H and P cases. In the homeotropic case 
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FIG. 1. (Color online) Definition of the homeotropic (a) and 
planar (b) geometries. They consist in a one-dimensional mul- 
tilayered structure made of Si02 (yellow slabs) and Ti02 (ma- 
genta slabs) in which is embedded a nematic liquid crystal 
layer. The uniform spatial distribution of its optical axis at 
rest is no = e z in the case H and no = e x in the case P. The 
excitation light beam impinges at normal incidence onto the 
structures, ko = fcoe z , and the incident polarization state is 
set to circular in the case H, Ei Ilc = E a (e x + ie y )/y/2, and to 
linear in the case P, Ei nc = Eoe y . 

with n c ff = n\\n±/(n 2 cos 2 0h + n\ sin 2 ©h) 1 ^ 2 , ^o = 
2tt/\q and Ao is incident wavelength in free space, 
whereas the condition 0p = 7r/2 in the planar case im- 
poses a constant phase delay 

A P = k Q L(n\\ - n±) . (7) 

Finally, we define the reduced incident light intensities 
as pu = /i nc / In in the homeotropic case, where /j nc is the 
incident intensity and /H=27r 2 cA3n 2 /[(n 2 — n^n^L 2 ] is 
the Freedericksz threshold (i.e. the intensity above which 
the initial director state no is unstable) for circularly po- 
larized light for a slab alone [16], whereas pp — I luc /Ip 
with 7p=87r 2 cA 2 ri|| (riii — ?i^)/[A 2 (?i|| +n±)] in the planar 
case [19], where c is the speed of light in free space. 

B. Light-matter interaction geometries 

A general feature of the optical reordering of liquid 
crystals is its strong dependence on the light-matter in- 
teraction geometry, namely the director field at rest and 
the incidence angle, polarization and intensity of the in- 
cident light field. Therefore a preliminary analysis is 
performed for the H and P cases from the knowledge 
of the reorientation dynamics for a nematic slab alone 
(i.e., without periodic structure) in order to identify sit- 
uations that will potentially exhibit a dynamical optical 
response. 

Homeotropic case - - The optical reordering of a 
homeotropic nematic film has been intensively studied 
for normal or oblique incidence and for a various set of 
incident beam polarization state, be it linear, elliptical, 
circular or unpolarized. Various static, periodic, quasi- 
periodic and aperiodic reorientation dynamics have been 



A h (t) = k a L [ [n eff (£, t) - n x ] d£ , (6) 
Jo 
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predicted and observed (sec [8] for a review). The cir- 
cular polarization case at normal incidence is a repre- 
sentative example for which periodic, quasi-periodic and 
chaotic rotations (see [20] and references therein) have 
all been reported. Therefore, it seems reasonable to ex- 
pect dynamical behavior too when a periodic structure is 
at work. However, we note that a thin nematic slab can 
be a crippling drawback since the light-induced dynami- 
cal richness declines as the thickness decreases, as shown 
in Fig. 2. The phase delay Ah, and the collective direc- 
tor precession angular velocity, Q = /dr, are plot- 
ted in Fig. 2(a) and 2(b), respectively, for two different 
thicknesses. From this figure, we see that the precession- 
nutation regime (i.e., f2 ^ with dApi/dr ^ 0) dis- 
appears below a typical cell thickness that corresponds 
to L = 2 — 3 /im for typical wavelength Ao = 600 nm. 
The typical reorientation dynamics diagram observed for 
thick cells (L ~ 100 /im, sec [17]) is thus qualitatively 
preserved in our case, where L = 5 /im. 

Planar case — The planar alignment geometry for ne- 
matics at normal incidence has been much less stud- 
ied than its homeotropic counterpart. To our knowl- 
edge, this geometry has only been considered in [21] 
and no further studies have been performed on it dur- 
ing the past two decades. This might be explained 
by a (twist) Freedericksz transition threshold predicted 
to be up to three orders of magnitude larger than the 
homeotropic case, thus preventing from experimental ob- 
servation. Moreover, the above-threshold director distor- 
tions have been predicted to be static whatever the light 
intensity [21], hence limiting the need for a deeper analy- 
sis. Nevertheless, very recently, the planar geometry has 
attracted a renewal of interest following seemingly exper- 
imentally accessible intensity threshold using a periodic 
structure [12]. 

In fact, the case of a planar slab alone was revisited 
in [19] where it is rigorously demonstrated that the reori- 
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FIG. 2. (Color online) H slab alone. Total phase delay Ah 
(a) and director precession angular velocity Q. — d<&^p /dr (b) 
vs. intensity for two different nematic thicknesses L = p\ 
(p — 5 and 10), where A = Xo/n± is the wavelength in the 
unperturbed nematic. Solid (dashed) curves refer to stable 
(unstable) states. The gray region in panel (a) corresponds to 
the range of values explored during the quasi-periodic regime. 



0.5 



S 3 o 



-0.5 




4 

R, 



FIG. 3. P slab alone. Output reduced third Stokes parameter 
S3 vs. intensity for three different nematic thicknesses that 
correspond to Ap/iv = 0.5 (a), 0.8 (b) and 1.6 (c). Solid 
(dashed) curves refer to stable (unstable) states. The gray 
regions in panels (b) and (c) correspond to the range of values 
explored during the oscillatory regime. 



entation threshold is significantly lower than previously 
predicted [21]. Moreover, a static distorted state is ex- 
pected above the reorientation threshold only for opti- 
cally thin enough liquid crystal layers, Ap/n < 0.64 [19], 
whereas a dynamical regime that takes place via a Hopf 
bifurcation is predicted otherwise, which was missed so 
far [21]. A secondary instability (heteroclinic bifurca- 
tion) has also been predicted, which eventually leads to 
a static regime, is associated (0.64 < Ap/n < 1), or not 
(Ap/n > 1), to a hysteresis behavior [19]. These re- 
sults are summarized in Fig. 3 where the reduced third 
component of the Stokes parameter at the output is 
plotted as a function of the reduced input intensity pp 
for three representative optical thicknesses that corre- 
spond to Ap/n = 0.5 [panel(a)], 0.8 [panel(b)] and 1.6 
[panel(c)]. 

Since L = 5 /im corresponds to Ap /n ~ 4 in the visi- 
ble range with the present choice of material parameters, 
a non trivial sequence of nonlinear reorientation dynam- 
ics is a priori expected. Importantly, note that the cho- 
sen film thickness is nevertheless thin enough to prevent 
from the influence of stimulated scattering, in contrast to 
thick cells where full polarization conversion can be ob- 
served [22, 23]. Indeed, the latter phenomenon is all the 
more important as the ratio A = L(nu — n±)/X between 
the thickness and the intensity modulation grating period 
arising from the coherent superposition of extraordinary 
and ordinary waves is large compared to unity. Since 
A ~ 2 in our case, the stimulated scattering can there- 
fore be safely neglected in our approach. 

To summarize, we further consider the case H under 
circularly polarized excitation at normal incidence and 
the case P under ordinary linearly polarized excitation 
at normal incidence, as sketched in Figs. 1(a) and 1(b), 
respectively. 
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III. SPECTRALLY MODULATED 
OPTOMECHANICAL EFFICIENCY 

In both H and P cases, the incident beam polariza- 
tion corresponds to an ordinary wave for the nematic 
at rest and the optical torque density exerted on no is 
thus zero. However the orientational state n = no is 
linearly unstable with respect to fluctuations when the 
light intensity is high enough to overcome the restoring 
elastic torque density that originates from the cell walls, 
which defines the optical Frecdericksz transition thresh- 
old [16]. The presence of the periodic structure does not 
suppress the optical destabilization scenario, but enriches 
it. Indeed, the reduced thresholds pH,th [Fig- 4(a)] and 
pp,th [Fig. 4(b)] strongly depends on the excitation beam 
wavelength Ao inside the bandgap, as already discussed in 
details in the homeotropic case [11, 13]. Namely, the de- 
fect modes evidenced by sharp peaks in the transmission 
spectrum Th [Fig. 4(c)] are associated to an enhanced 
light confinement in the nematic defect layer that lead to 
a significant reduction (up to several orders of the mag- 
nitude) of the required incident intensity to trigger the 
liquid crystal reorientation [Fig. 4(a)]. 

In the planar strongly modulated reorienta- 

tion threshold is observed too, as shown in [Fig. 4(b)], 
however, its fine structure is more complex than its 
homeotropic counterpart. Indeed the spectrum of /fp,th 
exhibits two sets of resonances with either symmetric [cir- 
cle symbols in Fig. 4(b)] or asymmetric [square symbols 
in Fig. 4(b)] lineshape that are related to the o- and e- 
defect mode wavelengths \° A and AJj, respectively. Such a 
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FIG. 4. (Color online) Reduced thresholds pH.th (a) and pp,th 
(b) and transmission spectra Th (c) and Tp (d) at rest (n = 
no). In panel (b), circle (square) symbols refer to the ordinary 
(extraordinary) mode shown in panel (d). 



behavior is a direct manifestation of competing influence 
of the o- and e-modes. On the one hand, symmetric reso- 
nance lineshapes in the planar case are reminiscent of the 
incident ordinary polarized light field enhancement in the 
nematic defect layer when Ao ~ A^, as in the homeotropic 
case. On the other hand, the origin of the asymmetric 
resonance lineshapes when Ao ~ A| can be qualitatively 
grasped by considering a fluctuation around the initial 
state n for a planar slab alone, see Fig. 5(a). The elas- 
tic torque density (r e i) driven by the cell walls tends to 
restore the director perturbation whereas the subsequent 
polarization changes (see typical polarization state inside 
the liquid crystal labeled as Elc in Fig. 5) are associ- 
ated to a spin angular momentum deposition [24] , hence 
an optical torque density (r opt ), that tends to increase 
the initial amplitude of the fluctuation, as illustrated in 
Fig. 5(a). The influence of the periodic structure is ad- 
dressed by noting that (i) the phase of the o-wave is al- 
most constant in the vicinity of a e-defect mode (except 
when A^j ~ X^) and (ii) a tt phase shift is experienced by 
the e-wave when incident beam wavelength explores the 
resonance linewidth. Therefore, the sign of the optical 
feedback (<5r opt ) is either negative or positive feedback 
depending on the resonance side, as shown in Fig. 5(b) 
and 5(c). We note here the analogy with the Fano reso- 
nance arising from the interference of resonant and non- 
resonant scattering waves, which is characterized by an 
asymmetric lineshape [25]. Indeed, resonantly excited e- 
wave leads to constructive and destructive interference 
with o-wave in the vicinity of a e-defect mode. When 
AJ] w AJ^j [i.e., A ~ 673 nm, see arrow in Fig. 4(d)], the 
dephasing associated to the defect modes between the 
o- and e- waves is even smaller than A^ is close to AJj, 
which eventually leads to symmetric resonance lineshape 
for pp.th [ see arrow in Fig. 4(b)]. 

The nature of the optical reordering transition also 
depends on the excitation beam wavelength. In the 
homeotropic case, the order of the orientational instabil- 
ity can be either first- or second-order depending on the 
detuning Aq — AJ] between the pump wavelength and the 
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FIG. 5. (Color online) Qualitative interpretation of the asym- 
metric resonance lineshape nearby an extraordinary defect 
mode in the planar geometry. Elastic and optical torque den- 
sities exerted on the director in presence of an orientational 
fluctuation are shown in panel (a) for a slab alone whereas 
panels (b) and (c) illustrate the negative and positive feed- 
back from the periodic structure for Ao < and Ao > AJ, 
respectively. 



5 



Hopf bifurcation 




Negative detuning 



Positive detuning 



640 650 660 670 680 690 700 710 

A. (nm) 

FIG. 6. Spectrally modulated nature of the reorientation 
transition that corresponds to a stationary bifurcation except 
in the vicinity of e-defect modes, which is identified by gray 
regions, where a Hopf bifurcation takes place. 



nearest defect mode A^, as previously discussed for iden- 
tical structure with linearly polarized excitation [11, 13], 
which is a general result for homeotropic films inside a 
resonator under linearly polarized light [26, 27]. Indeed, 
the transition is first-order for positive detuning whereas 
it is second-order for negative detuning [see later Fig. 8(a) 
and 8(b) in the next section]. Such a behavior is there- 
fore independent on the presence or not of an intrinsic 
optical bistability (i.e., for the nematic slab alone) since 
linearly and circularly polarized excitations are related 
to second- and first-order transitions, respectively. In 
the planar case, we also find that the order of the tran- 
sition is controlled by the detuning with respect to the 
nearest o-defect mode, whatever the detuning between 
the o- and e-defect modes, but with some exception in 
the vicinity of the e-defect modes where a supercritical 
Hopf bifurcation is found (i.e., second-order transition), 
as illustrated in Fig. 6. This emphasizes the prominent 
role played by the e-defect modes in the case P although 
incident light is ordinary polarized. 
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FIG. 7. (Color online) Director regimes in the H (a) and P 
(b) geometries. Uniform precession of the director around the 
z axis is found in the case H whereas either static distorted 
states (see markers) or director oscillation in the (x, y) plane 
are predicted in the case P. 
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FIG. 8. Director reorientation characteristics for negative 
[panels (a,b,c)] and positive [panels (d,e,f)] detuning Ao — A^ = 
±5 nm near the defect mode A^ « 600 nm. The total phase 
delay An/2n and angular rotation frequency Q vs. input in- 
tensity are shown in panels (a,d) and (b,e), respectively. The 
total twist defined by * = $h(1,t) — $h(0,t) vs. Au/2n 
is shown in panels (c,f). Solid (dashed) curves correspond to 
stable (unstable) regimes. 



IV. ALL-OPTICAL DYNAMICAL OPTICAL 
RESPONSE 

As expected from the analysis of H and P nematic 
slab alone in section II B, we find dynamical regimes for 
the director reorientation above /?H,th and yOp,th, respec- 
tively. In the case H, however, there is no nutation regime 
(<96h/<9t 7^ 0) whatever is the wavelength although uni- 
form precession (d^n/dr = constant) is always present. 
In the case P, the distorted state is static or correspond to 
an oscillatory motion in the (x, y) plane. These regimes 
are sketched in Fig. 7 and next we discuss the optical 
response of the structure when dynamics takes place. 

Homeotropic case — The uniform director nutation- 
free precession dynamics corresponds to a constant polar 
spatial profile (Oh) for which Ah is an estimate of the 
overall amplitude [see Eq. (6)]. The azimuthal spatial 
profile is also constant but in a frame that rotates at con- 
stant angular velocity VL = d^n/dr = d&^p /dr around 
the z axis. Also, we define the overall twisted charac- 



tcr of the distortion profile in the rotating frame as the 
twist angle between the output and input facets of the 
nematic defect layer, '4? = <I>h(1,t) — 3>h(0, t). The in- 
cident intensity dependence of these three representative 
quantities, Ah, ^ and *S? are shown in Fig. 7 for negative 
and positive detuning, Ao — A^ = ±5nm, near the defect 
mode A^ re 600 nm [see Fig. 4(c)]. 

In fact, Ah and fi are related to distinct features of the 
polarization state at the output of the structure, namely 
the ellipticity and the rotation rate of the polarization el- 
lipse, whose azimuth is given by <f>n(l,T) (non-adiabatic 
light propagation when '4? ^ slightly modify this corre- 
spondence) [17, 24]. Since Ah is time-independent, the 
total transmission Th is constant and the dynamical op- 
tical response of the structure is merely contained in the 
polarization azimuth whatever the sign of the detuning. 
The intensity dependence of the optical response, how- 
ever, strongly depends on wavelength, as shown in Fig. 8. 
Indeed, for negative detuning, the overall reorientation 
picture is very similar to the case of linearly polarized 
excitation [11], as shown in Fig. 8(a). This is due to 

the almost untwisted distortion profile, ^h"" 1 -C 1, hence 
$ < 1, see Fig. 8(c). This is in stark distinction with 
the case of positive detuning where a complex reorienta- 
tion diagram is found [Fig. 8(c)] as a result of large twist 
amplitude and associated significant non adiabatic light 
propagation effects. Such a behavior is not observed for a 
homeotropic slab alone and is the signature of the inter- 
action between o- and e-defect modes that are no longer 
degenerate when no. 

Planar case — Static or oscillatory regimes are found 
depending on the incident wavelength and intensity. 
Static distortions are found for almost any wavelength ex- 
cept in the vicinity of extraordinary defect modes where 
oscillatory dynamics appears via a Hopf bifurcation, as 
illustrated in Fig. 9(a). In that case, the output optical 
response exhibits time-dependent o- and e-transmissions, 
as shown in Fig. 9(b), however noting that the modula- 
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FIG. 9. (a) Reduced Stokes parameter S3 vs. input intensity 
where solid (dashed) curves refer to stable (unstable) states, 
(b) Ordinary and extraordinary transmissions T (solid curve) 
and T e (dashed curve) vs. input intensity at Ao = 650nm. 
Gray areas correspond to the range of values explored during 
the dynamics. 
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FIG. 10. (Color online) Fourier spectra of the time-dependent 
total transmission [panels (a) and (b)] and spatio-temporal 
dynamics of the azimuthal spatial profile <J?p(£, r) [panels (c) 
and (d)] just above the Hopf bifurcation threshold [see A in 
Fig. 9(a)] and far above the threshold [see B in Fig. 9(a)]. 



tion depth of T (r) are much less pronounced than the 
one of T e (r) and T e <C T . Details of the dynamics is il- 
lustrated in Fig. 10 near and well above the threshold at 
intensities that correspond to labels A and B in Fig. 9(a), 
respectively. As expected from a Hopf bifurcation, a sin- 
gle frequency characterizes the dynamics near pp.th, see 
Fig. 10(a), whereas a second frequency grows up as the 
intensity increases, which leads to a quasi-periodic oscil- 
lation with the subsequent appearance of harmonic fre- 
quencies and frequency mixing that are clearly seen in the 
Fourier spectrum of the time-dependent transmission, see 
Fig. 10(b). The corresponding spatio-temporal director 
dynamics, <&p(£,r), is shown in Fig. 10(c,d). Aperiodic 
dynamics, which has not been predicted for a planar slab 
alone [19], is also found in simulations at higher intensi- 
ties (not shown here), which emphasizes how the presence 
of the periodic structure contrary to the homeotropic 
case. 



V. CONCLUSION 

The dynamical response of all-optical liquid crystal in- 
filtrated photonic structure has been discussed in two 
light-matter interaction geometries that correspond to a 
circularly [linearly] polarized incident beam impinging at 
normal incidence onto a one-dimensional periodic struc- 
ture in which is embedded a nematic liquid crystal defect 
layer with homeotropic [planar] alignment (i.e., perpen- 
dicular [parallel] to the layer) . In both cases the incident 
beam polarization state corresponds to an ordinary wave 
and optical reordering takes place only above a thresh- 



old for the input intensity that is much smaller than the 
usual optical Freedericksz transition threshold reduced 
with respect to in the vicinity of ordinary defect modes. 
It has been shown that self-sustained dynamics of the 
output light polarization and/or intensity can take place 
in contrast with all previous attempts discussed so far, 
thus paving the way towards all-optical dynamical re- 
sponse of photonic structures using liquid crystals. For 
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